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Abstract 

The renormalization group is not only a powerful method for de- 
scribing universal properties of phase transitions but it is also useful 
for evaluating non-universal thermodynamic properties beyond mean- 
field theory. In this contribution we concentrate on these latter aspects 
of the renormalization group approach. We introduce its main under- 
lying ideas in the familiar context of the ideal Bose gas and then apply 
them to the case of an interacting, confined Bose gas within the frame- 
work of the random phase approximation. We model confinement by 
periodic boundary conditions and demonstrate how confinement mod- 
ifies the flow equations of the renormalization group changing thus the 
thermodynamic properties of the gas. 



1 Introduction 



The recent experimental realization of Bose-Einstein condensation has re- 
newed the interest in the study of trapped, dilute, weakly interacting Bose 
gases. Despite the long history of this subject there still exist unresolved 
theoretical problems for the interacting Bose gas even without the presence 
of a trap. Here we discuss some of these problems, their origins, and possible 
ways of overcoming them. 

A trivial example of Bose-Einstein condensation occurs, of course, in the 
non-interacting (ideal) case. What is of real physical interest though is the 
interacting case. This is so not only because the ideal g 8bS IS, clS the name 
suggests, an idealized approximation to reality but also because it is the 
interactions that render the condensate a true superfluid with remarkable 
properties . The presence of interactions, even weak ones, complicates the 
mathematical treatment of the problem and calls for approximate schemes. 
The most well-known is the Bogoliubov approach and, though it was pro- 
posed more than 50 years ago, it remains the approach used most frequently 
in the literature of the subject. 

The first step of the Bogoliubov approach is to simplify the interparticle 
interactions. We consider that the origin of the interactions is two-body 
collisions only. The gas is assumed to be very dilute, which is equivalent to 
assuming that the interactions are very weak. This means that the interacting 
gas differs only slightly from the ideal gas and consequently most particles 
in the interacting ground state should have zero momentum, as they do 
in the ideal case. In this limit of low momenta, one can show that the 
interactions are dominantly of s-wave type. Thus they can be characterized 
by a scattering length a which is positive for repulsive interactions. 

The second step is to replace the creation and annihilation operators of 
particles in the zero- momentum state Og and oq, respectively, by their expec- 
tation values. The reason is that we will want to treat the interaction term 
perturbatively. However, for bosons, the usual form of perturbation theory 
fails, because neither Qq nor annihilate the interacting ground state. This 
means that it is impossible to define normal-ordered products with vanishing 
ground-state expectation value and consequently Wick's theorem, a corner- 
stone of conventional perturbation theory, cannot be applied in a straight- 
forward way This problem is fixed if one replaces Oq — > (oq ) = ^/No and 
ao —>■ (cto) = \/Nq, where Nq is the average number of condensate atoms, 
and reinterprets Oq", as ordinary numbers that commute. Approximating 
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these operators this way is justified because their commutator is of order 
1/\/Nq with respect to their individual matrix elements and because we as- 
sume that the condensate density uq = Nq/V, V being the volume of the 
system, remains finite as A^^o, — > oo in the thermodynamic limit. In fact, 
this replacement leads to the so-called, Hartree-Fock-Bogoliubov (HFB) the- 
ory 1^. One can do better than this by replacing the operators by their 
average value (which is an ordinary number) plus a small fluctuation (which 
is an operator), i.e. Oq — >■ \/No + Sqq, Qq + 6ao, and minimizing the 

Hamiltonian with respect to the number of condensate atoms Nq. This leads 
to the Bogoliubov theory proper (B) 0, whose results, of course, reduce to 
those of HFB when the fluctuations are set to zero From these con- 

siderations it is apparent that the Bogoliubov theory is of mean-fleld type, 
so one could in principle improve upon it by using more sophisticated tech- 
niques. One promising possibility, near the critical region, is renormalization 
group @, 

What makes the use of beyond mean-field theory methods imperative 
though is that the Bogoliubov theory has other shortcomings in addition to its 
mean-field character 0. In the critical region the Bogoliubov theory simply 
does not work because there are fiuctuations around the mean-field that 
cannot be treated perturbatively. This happens because, as the temperature 
approaches the critical temperature Tc, the thermal cloud density develops an 
infrared singularity and thus diverges as the momentum tends to zero p, |[ . 
This really calls for Wilsonian renormalization treatment [|TU| of the problem 



around the critical region. It is well-known that renormalization techniques 
circumvent such infrared singularities because, to put it simply, one arrives 
at the fiow equations without having to integrate out all momenta down to 
zero 
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Furthermore, within the Bogoliubov theory, it is usual to make a further 
approximation. One assumes that the presence of the condensate and of the 
thermal cloud does not modify significantly the effective interaction between 
two colliding atoms from its vacuum value. This is the Bogoliubov-Popov 
theory (BP) [0. To avoid this approximation, which is particularly unjus- 
tified in two-dimensional cases or when the interactions are attractive, one 
has to use many-body T-matrix theory [^, However, this approach 

also runs into infrared singularity problems in the critical region. Therefore 
renormalization techniques would have to be employed to take correctly into 
account the many-body effect on two-body collisions near the critical region. 

However, the renormalization group is not only a powerful method for 
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describing interacting quantum systems close to a phase transitions but it 
is also useful for evaluating thermodynamic properties at arbitrary temper- 
atures thereby transcending mean-field theory. In this contribution we con- 
centrate on these latter aspects of the renormalization group approach by 
introducing its main underlying ideas in the familiar context of the ideal 
Bose gas and by applying them to the case of an interacting, confined Bose 
gas within the framework of the random phase approximation. 

This paper is organized as follows: In section 2, we present an introduction 
to renormalization group methods in the familiar context of the ideal Bose 
gas, in an arbitrary number of spatial dimensions. In section 3 we apply the 
renormalization group to the realistic case of an interacting, confined Bose 
gas whose confinement can be modelled by periodic boundary conditions. 
Comparisons with the ideal gas are made where appropriate. We analyse 
the resulting critical fixed point and its associated critical exponent char- 
acterizing the scaling of the correlation length in three spatial dimensions. 
We also investigate a particular non-universal property which is accessible 
to experimental observation, namely the second order coherence factor, and 
discuss the effect of one-dimensional confinement on this physical quantity. 

2 Renormalization Group and the Ideal Bose 
Gas 

Before treating the interacting Bose gas case, we illustrate the basic ideas of 
renormalization group in the familiar context of the ideal Bose gas. 

Let us start from the path integral representation of the grand-canonical 
partition function of the ideal Bose gas in D dimensions [^, i.e. 

with the (dimensionless) action 

SU), 0*1 = - / dT rf^x 
n Jo Jv 



d h 



(^'^)[^^-^V2-/x]0(r,x) 



The boson mass is denoted by m, /i is the chemical potential and P = l/^ksT) 
is the inverse temperature with Boltzmann's constant ks- Due to the trace 
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operation involved in the definition of tlie partition function the complex- 
valued field 0(r, x) has to fulfill periodic boundary conditions with respect to 
the imaginary time r, i.e. 0(r, x) = 0(r + /i/3,x). Thermodynamic properties 
of the non-interacting Bose gas can be evaluated from the extensive quantity 



f3,V) = - In /3, V) = pn{fx, P,V) =Y^ In 



1 _ g-z^C^-^) 



(2) 



with V denoting the volume of the system. f3,V) is the grand-canonical 
thermodynamic potential. Modelling confinement by periodic boundary con- 
ditions with spatial periods Li (i G {1,...,Z)}) the corresponding possible 
wave numbers are given by kj = 27rnj/Lj with rii being integer. 

If we want to evaluate the thermodynamic function of Eq.(Q) numeri- 
cally, it is convenient to introduce first of all an ultraviolet momentum cutoff 
l^kmaxl = which constitutes an upper limit for the summation over all 
states of the ideal Bose gas. In order to ensure convergence of the partition 
function this momentum cutoff has to be much larger than all physical mo- 
menta characterizing the problem, i.e. (/iA)^/(2m) S> 1/(3, ii. Introducing 
this ultraviolet momentum cutoff and the dimensionless scale parameter / 
according to the relation |k| = Ae~', in the continuum limit, i.e. Lj oo, 
Eq.(^ can be approximated by 

F(Mo, bo, V) ^ r dl' M e-^'' In \l - e-^(^'){|-M(/'))l (3) 

with the dimensionless, rescaled thermodynamic parameters 

h{l) = 60 ^ (/5//5a) 
M{1) = Moe^' = iMe^' (4) 

and with the cutoff dependent scale factor P\ = m/ {TiKY. The quan- 
tity Mdl' = [VA^Qo/i'^'^)^] dl' denotes the number of states in volume 
V within the infinitesimal momentum shell ^A(l — dl') < |^k| < ^A. = 
271^^'^ /T(D /2) is the surface of a D- dimensional hypershpere of unit radius 



and r(x) denotes the Gamma function |TH . 

The main idea underlying the renormalization group approach is to per- 
form the integration over all states involved in the evaluation of the partition 
function or of the thermodynamic function of Eq.(|^) not in one step but in 
small, successive steps. Integrating out some of the states in one of these 
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small steps yields scaling relations from which one can determine the non- 
analytic behaviour of thermodynamic functions at a phase transition. Suc- 
cessive applications of this partial decimation procedure eventually yields the 
partition function. 

In order to derive such a scaling relation for the thermodynamic function 
of Eq.(|^) let us integrate out a small momentum shell corresponding to the 
interval /' G [0, /], for example. For the ideal Bose gas this yields the result 

F(Mo, 60, V) = {f dl' + r dl') {M e-™' In [l - e"''('')(5-^^(n)l } 

Jo Jl L J 

dl' M e-™' In fl - e^^('')(^-*^(''»l + e"^' F(M(/), V). (5) 





The last equality of Eq.(|]) has been obtained by shifting the integration vari- 
able /' in the second integral according to — /. The two characteristic 
steps of Wilsonian renormalization are apparent in the derivation of this scal- 
ing relation. The first step, the Kadanoff transformation, involves a partial 
decimation of some states of the Bose gas. In the second step the original 
ultraviolet momentum cutoff UK is re-established by rescaling momenta ac- 
cording to the transformation |^k| |^k| eK In our case, this latter trivial 
rescaling is performed by the translation /'—;>/' — /. In configuration space 
this trivial rescaling implies that we are increasing the minimum distance of 
resolution or effective block size from 1/A to 1/(A e^') so that with increasing 
values of / the physical system is described on larger and larger length scales. 
Thus, expressed in terms of this minimum distance of resolution any physical 
length L shrinks with increasing values of / according to L — Le~K Scaling 
relations of the form of Eq.@ are useful for investigating the non-analytic 
scaling behaviour of thermodynamic quantities close to a second order phase 
transition p. 

Applying the scaling relation of Eq.@ repeatedly the evaluation of the 
thermodynamic function F{Mq, bo, V) can be reduced to the solution of a 
system of ordinary differential equations, i.e. 



g-6(Z)(i-M(0) 



dl 
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which have to be solved with the initial conditions M{1 = 0) = Mq <^ 1 
and b{l = 0) = Bq ^ 1 consistent with the choice of the ultraviolet cut- 
off. Thereby, the thermodynamic function F{Mo, bo, V) = F{1 oo) is 
obtained from the solution of Eqs.(|^) with the additional initial condition 
F{1 = 0) = 0. We will be referring to Eqs.(^ as the renormalization group 
(RG) equations. We note that the equation for F{1) depends on the solutions 
M{1) and h{l) but F{1) itself does not couple back to the equations of motion 
for these two quantities. Thus, one can first solve the two equations of mo- 
tion for the scaled thermodynamic parameters M{1) and h{l) and insert the 
resulting solutions into the equation of motion for F{1). In the case of the 
ideal Bose gas the RG equations for M{1) and h{l) are simple and describe 
the trivial scaling of these scaled thermodynamic parameters. This example 
demonstrates the basic concepts involved in the evaluation of thermodynamic 
partition functions by RG methods. 

The system of Eqs.(^ is autonomous with an unstable, fixed point at 
(6*,M^,) = (0,0). The parametric plot of a typical trajectory M{1) versus 
h{l) is drawn in Fig.^ together with this fixed point and its associated stable 
(s) and unstable (u) manifolds. It is well known that in the grand-canonical 
ensemble the phase transition of the ideal Bose gas occurs at zero chemical 
potential, i.e. at Mq = 0. According to Eqs.(|D this reflects the fact that 
at the critical point thermodynamic properties of an ideal Bose gas are de- 
termined by the stable manifold of the unstable fixed point (6*,M*). The 
unstable manifold of this fixed point governs the behaviour of the ideal Bose 
gas close to criticality. The eigenvalues of the scaled thermodynamic parame- 
ters corresponding to the stable and unstable manifolds are given by A_ = —2 
and A+ = 2 (compare with Eqs.(^)). In particular, the positive eigenvalue 
A_|_ = 2 which is associated with the (one dimensional) unstable manifold 
determines the critical exponents governing the scaling relations of singular 



parts of physical quantities near the critical point |rT|. As an example, 
let us consider the correlation length of the ideal Bose gas which changes 
under renormalization according to ^(/) = Furthermore, according to 

Eqs.(§), the relevant variable scales as (M(/) - M*) = (Mq - M^)e^+K If the 
system is close to the critical point, i.e. | Mq — l-C 1, and if we iterate 
the renormalization group equations up to a point ^ 1 far away from this 
critical point, say with M{Iq) — = 1, we find the scaling relation 

eo = = miMo - M,)-i/"+ = e(/o)(Mo - M,)-^ (7) 

Thus, the divergence of the correlation length at the critical point is governed 
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by the critical exponent v = 1/A+. With the help of scaling relations such as 
Eq.(^, one may also relate the characteristic eigenvalue A+ of the unstable 
manifold to the critical exponents of other physical quantities of interest. Let 
us consider the thermodynamic function F{Mq, bo, V) as a further example. 
Whereas for any finite value /q the integral in the last line of Eq.(^ is a 
smooth function of the intensive thermodynamic parameters fi and (3, the 
second term of this last line gives rise to a non-analytical behaviour in the 
neighbourhood of the critical point. If we are close to the critical point and 
choose Iq again so large that (M(/o) — M^,) = 1 and 6(/o) -C 1 we find that 
the singular part of this thermodynamic function scales as 

Fs{Mo, bo, V) = e-™° 0, V) = (Mo - M,)^/^+ 0, V). (8) 



3 The Interacting Bose Gas 

Treating two-body collisions between bosons in the low-momentum or s-wave 
approximation the path integral representation of the partition function of 
the homogeneous interacting Bose-gas is given by 

/3, V, g) ^ Tre-^(^-^^) = j 5[<P, rje'^l*'^*] (9) 

with the (dimensionless) action 



S[(P, 0*] = - / dr d^^ 



h Jo Jv 



(10) 

In the low-momentum approximation the interparticle interaction can be de- 
scribed by the zero-momentum component of the Fourier transform of the 
two-body interaction potential. Thus, within this approximation a repulsive, 
short-range potential can be characterized by a positive interaction strength 
g. In three spatial dimensions, for example, this interaction strength is re- 
lated to the positive scattering length a of the interparticle interaction by 
the familiar relation g = Anh a/m. 

In the rest of this paper we will be interested not only in the case of an un- 
confined interacting Bose gas but also in cases where at least one of the spatial 
degrees of freedom is confined by a potential. In the simplest approximation 
such a confinement in the 3- or 2;-direction, for example, can be described by 



7 



a periodic boundary condition of the form 



with Lz denoting the characteristic length of confinement. Within such a de- 
scription of confinement the translational invariant character of the problem 
is conserved. 

Before addressing the evaluation of the partition function of Eq. (|^) with 
the help of the renormalization group let us briefly summarize its approx- 
imate evaluation within the framework of mean-field theory. In mean-field 
theory one approximates the partition function of Eq.(^ by saddle point in- 
tegration []T^ . For this purpose one determines first of all the most probable 
configuration 0(x) by minimizing the action of Eq . (p!oD . This way one arrives 
at the Gross-Pitaevski equation M 



(x) = 0. 



(11) 



In the homogeneous case that we are examining, the most probable static 
and space-independent configuration is given by 



(12) 



In a second step one expands the action of Eq.(|TD|) around this most probable 
configuration up to second order assuming that the fluctuations around this 
most probable configuration are small. The resulting Gaussian integrations 
can be easily performed, thus yielding the well-known mean-field approxima- 
tion for the partition function [|18| . 



3.1 Renormalization group approach 

Using the renormalization group it is possible to improve on this mean-field 
approximation. However, contrary to the case of the ideal Bose gas, in the 
presence of interparticle interactions the first step of the renormalization 
group method, namely the Kadanoff transformation, can be implemented 
only approximately. A frequently employed approximation in this context is 
the random phase approximation |19| whose resulting renormalization group 
equations will be discussed in the following. 
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Kadanoff Transformation 



As we want to take into account tlie fluctuations of the complex-valued fleld 
(f){T, x) beyond mean-field theory, it is convenient to separate this field ac- 
cording to 

0(r,x)-.0 + </)(r,x). (13) 

Thus, we obtain the following symmetry broken form of the (dimensionless) 
action 

S[<P, 0*] = -pVifxno -^]+T dx 0*(a;) [/^ — - — - /i + 2gn,]<P{x) 



2 h J dr 2m 

2h 



I dx [<l>*{x)<f>*{x) + <P{xmx)] 



J [0*(a;)0*(x)0(x) + 0*(x)0(x)0(x)] 
+ J- J dx (f)*{x)(j)*{x)(j){x)(j){x) (14) 

with the condensate density uq = |0p = fi/g. The U(l) global symmetry 
of the action of Eq.([TO|) has been broken spontaneously by the introduction 
of the most probable configuration (p. In order to simplify the notation 
we have taken x to represent both time and coordinates, so, for example, 
(f){x) = 0(r, x) and dx = dr d^x. 

Now we want to integrate out large momentum fiuctuations, so we also 
split the fiuctuation 0(a;) around the most probable configuration into a 
long wave length component 0< {x) and into a short wave length component 
6(f>y{x), i.e. 

0(x) = 0<(x) + (50>(a;) (15) 



with 



and 



(I><{X) = E E -iTT "Pnm (16) 



where '{>nxa (^V'nm) are the Fourier components of 0<(a;) ((50>(a;)), no; = 
n {2n l(5K) are the Matsubara frequencies, and TiVk denotes a hypersphere of 
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radius hK in D-dimensional momentum space. The short wave length fluctua- 
tions involve momentum components which are contained only in an infinites- 
imally thin shell in momentum space of thickness ^A(l — dl) < \Kk\ < hA 
and which is denoted h6Vk (compare with Fig.|^). For the sake of simplicity 
we will be referring to (p<{x) as the lower field and to 50>(x) as the upper 
field. 

The upper field which we want to eliminate from the partition function 
involves an infinitesimal momentum shell. Therefore, it is sufficient to ex- 
pand the action of Eq.(0) up to second order in terms of the upper field 
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The resulting expression can be further simplified in the random phase ap- 
proximation [IS]. In this approximation it is assumed that the upper field is 
rapidly varying in space in comparison with the lower field so that the dom- 
inant contributions to the action arise from those particular terms quadratic 
in the upper field which are slowly varying in space. Thus, in the random 
phase approximation, the action reads 



^[0,0* 



S'[0<,0<] + - I dx 5(f)*y{x) 



h 



Or 2m 



-/i 



(50>(x) 



+ 



2n 



+ 



9_ 

2h 



dx |0 + 0<(a;)p 



dx (0 + 0<(x))" 



dx (0 + 0<(x))^ 



Vh(3 
1 

1 



dx 50^(a;) 50>(a;) 
dx 6(f)y{x)S(j)y{x) 
dx 6(f)y{x)6(j)^{x) 



:i8) 



where we have taken into account that terms linear in (50>(x) and 6(j)Z,{x) 
vanish because the fluctuations are expanded around the most probable con- 
figuration which fulfills the Gross-Pitaevski equation ([TTD . Performing the 
Gaussian integration in terms of the upper fields yields an effective action 
for the lower field of the form IHI 



'eff 



0*<] = ^[0<,0*<] + 5^[0<,0: 



(19) 



We want to point out that, contrary to the case of the mean-field approx- 
imation, now the contribution to the effective action, i.e. 55'[0<,0^], still 
depends on the lower field which characterizes the long wave length fluctu- 
ations. Stated differently, the mean-field approximation would correspond 
to the replacement 5S'[0<,0^] <55'[0< = 0,0^ = 0] in this step of the 
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decimation procedure. Thus, after integration over the short wave length 
fluctuations the simplest improvement transcending the mean-field approxi- 
mation is obtained by expanding 6S[(j)^, (p*^] up to second order in the lower 
field. As apparent from Eq.([T^, such a second order expansion leads to a 
change of the most probable configuration, i.e. 0^0 + 6(j), through the 
terms linear in 0< and (p*^ and to a change of the chemical potential through 
the terms quadratic in the lower field, i.e. fi —>■ 6fi. Enforcing the relation 
([T^ ) at each step of the renormalization procedure implies a correspond- 
ing scaling of the interparticle coupling strength g according to the relation 

Trivial rescaling 

The second step of the renormalization procedure is to recast the effective 
action of Eq.(|TUp in terms of the new chemical potential /ig// = fi+Sfi and the 



new interparticle coupling strength (/g// = 9 + ^9 the form of the original 
action by re-establishing the original momentum cutoff ^A. For this purpose 
we have to rescale momenta according to |k| — > |k(Z)| = \k\eK As in the 
ideal g cLS CcLSG, this transformation and the demand that the action remains 
formally the same after each renormalization step induce the appearance of 
rescaled parameters in the action, i.e. 



V - 


- vii) = 


re-™ 


p - 


- m = 




4> - 


- 0(0 = 




fi - 


-> M0 = 




9 - 


- 9(1) = 


(9 + Sg) e(2-^)' 



(20) 

These scaling relations reveal two differences from the corresponding scaling 
relations of the ideal gas. The first is that there is, of course, an additional 
equation for the rescaling of the coupling constant. The second is that the 
rescaled chemical potential contains an additional correction 6fi which itself 
depends on / and is therefore a non-trivial contribution to the scaling of 
This non-trivial part was not present in the ideal g We also note that 

the coupling constant equation contains a non-trivial part 5g as well. 
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The Confined Renormalization Group Equations 



Because we have taken the renormahzation step dl to be infinitesimal, the 
corrections for the chemical potential and the coupling constant take the 
form of differential equations with respect to the continuous parameter I. 
Introducing the dimensionless, scaled coupling constant 



G{1) = (3a g{l) 



(21) 



and the scaled thermodynamic parameters of Eqs.(§) we obtain the renor- 
malization group equations for the interacting Bose gas M 



dF{l) 

dl 
dM{l) 

dl 
dG{l) 

dl 

dbjl) 
IT 



VA^ d{l) e 



-Dl 



In (2sinh[A(0/2]) - 6(/)(e> - Mq e'')/2 
2M{l) + d{l) G{1) A[M{l),h{l)i 
-{D - 2)G{1) + d{l) [G{l)f B[M{1), 
-26(/) 



(22) 



with 



A[M{l),b{l)] = b{l) [2M(/) -2e>] - 



2A(/) 
1 



B[M{l),b{l)] = 3b{l) 



2sinh2(A(/)/2) 
coth[A(/)/2] 



+ 



2A(/) 



coth[A(/)/2] 

- im? ^ 



M(/) 
2A(02 

[2e> + M(/)]2, 



+ 



2A(/)2 
coth[A(0/2] 



_2sinh2[A(/)/2] ' A(/) 
and with the scaled cutoff energy e> = 1/2. The quantity 

X{l)=b{l)^e^ [6> + 2M(/)] 



[2e> + M{l)Y (23) 



(24) 



indicates that at each step of the renormalization procedure energy and mo- 
mentum are related by the Bogoliubov dispersion relation. 

Contrary to the corresponding renormalization group equations of the 
ideal gas (compare with Eqs.(^) Eqs. (^2D involve two nonlinear coupled 



12 



equations of motion for the scaled parameters M{1) and G{1). The corre- 
sponding mean-field results would be obtained by setting the fluctuations 
around the mean-field equal to zero. In this case, both A[M{l),b{l)] and 
B[M{l),b{l)] vanish and the equations of motion for M{1) and G{1) decou- 
ple. Furthermore, comparison with Eq.(P) reveals that in the mean-field 
approximation the equation for the scaled chemical potential M{1) is identi- 
cal to the corresponding equation for the ideal gas which just describes the 
trivial scaling of M{1). It should be mentioned that, contrary to the trivial 
scaling of M{1), the trivial scaling of the dimensionless coupling parameter 
G{1) depends on the dimensionality of the problem. 

The quantity d{l) appearing in these RG equations describes effects of 
confinement as long as they can be modelled by periodic boundary conditions. 
These periodic boundary conditions lead to a quantization of the momentum 
components in the confined directions. Thus the free-space integrals over 
momentum which appear in the continuum limit are replaced by sums of the 
form 

J2 =V{Ae~Yd{l)dl + 0[{dl)% (25) 

where d{l) characterizes the number of momentum states which are con- 
tained in an infinitesimal momentum shell whose thickness is proportional to 
dl (compare with Fig.^). For example, in the case of confinement of one spa- 
tial direction, say the z-direction, we obtain in the case of a D-dimensional 
problem the expression 

^ ' (27r)^ L,e-^A ^ ' 

where denotes the length of confinement, ([x] denotes the largest integer 
which is less or equal to x.) Taking the limit — > cx), Eq. (p6[) yields 
d{l) = VLoI{'2.'k)^ and we obtain again the free-space renormalization group 
equations. For the special case of three spatial dimensions these free-space 
results reduce to those of Bijlsma and Stoof |^. 



3.2 Universal Critical Properties 

For an investigation of the universal critical properties of the interacting Bose 
gas we have to determine the unstable fixed point of the renormalization 
group equations Eqs. (^21) and its corresponding linearized flow equations. 
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For this purpose we have to study Eqs.(^) in the hmit Z — > oo in which the 
effective temperatures are large, i.e. /?(/) = Pq e~^' ^1 |^. In this hmit the 
quantities A[M{l),b{l)] and B[M{l),b{l)] simphfy to the expressions 

Alum h(lM ^ + ^^^^^'> + 2^'(0^> + 2M=^(06> 

A[M{l),b[l)\ -2 6(/)e|[e> + 2M(/)]2 ' 

mM(l)h(l)] 56| + 26>M(/)+2M^(0 

B[Mil)Ml)] b{l)ei[e,+2M{lW ' ^ ^ 

Therefore, the equation of motion for M(/) does not couple directly to G{1) 
but to the modified, scaled coupling G{1) = d{l)[G{l)/b{l)]. In the absence 
of any confinement, i.e. for d{l) = Q£)/[27r)^ , the trivial scaling of this 
scaled coupling is given by G{1) = G{0) e~*^^~^)' so that this quantity is 
relevant for dimensions D < A contrary to the variable G{1) itself whose 
trivial scaling indicates that it is relevant for dimension D < 2. However, in 
the case of confinement this trivial scaling changes. If one degree of freedom 
is confined, for example, Eq.(]26|) implies that for sufficiently large values of 
/ the trivial scaling of the scaled coupling is given by G{1) = G{0) e~(^~^~^)' 
which typically has a significant influence on the position of the unstable 
fixed point of the renormalization group equations. 

Let us consider the physically important case of D = 3 in the absence of 
confinement in more detail. It is straightforward to determine the unstable 
fixed point for the three dimensional case, i.e. 

= 1/2 

G, = {G/b), = 7rV2 (28) 

with G{1) = G{l)lb{l). The parametric plot of a typical trajectory G{1) versus 
M{1) is drawn in Fig.|^ together with the fixed point and the associated stable 
(s) and unstable (u) manifolds. The eigenvectors corresponding to the stable 
and unstable manifolds are obtained by linearizing the RG equations Eqs. (p2D 
around the fixed point. 



dl TT^ 



dAG{l) 27r2 
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with AM = M{1) - and AG = G{1) - G^. The eigenvalues and eigen- 
vectors of this set of hnear equations are given by 

A± = (3 ± V33)/6 (29) 

and 



(AM,AG')± = (^^,1). (30) 

According to Eq.(^ the positive eigenvalue A+ = (3 + -\/33)/6 describes the 
rate of increase of the relevant variable and determines the critical exponent 
for the correlation length, i.e. 

V = 1/A+ = ^ 0.686. (31) 

' 3 + V33 ^ ^ 



as already shown in . This critical exponent compares well with the known 
result of = 0.67 The corresponding mean-field value for the critical 
exponent is z/mf = 1/2, exactly the same as for the ideal gas. From these 
considerations it is apparent that taking into account fluctuations beyond 
the mean-field approach is crucial for describing the critical behaviour of the 
system. 

3.3 Non-Universal Critical Properties 

The renormalization group equations can also be used to calculate non- 
universal properties of the interacting Bose gas at the critical temperature, 
for example. For this purpose one has to solve the renormalization group 



equations of Eq. (|22D along the stable manifold of the unstable, fixed point 
(M*,^*). As an example let us consider the critical behaviour of the (spa- 
tially averaged) second order coherence factor g'^'^^ in the physically inter- 
esting case of D = 3 p]. This second order coherence factor is well-known 
in quantum optics from the correlation experiments of Hanburry-Brown and 
Twiss pB| on electromagnetic radiation from distant stars. It describes 
the (spatially averaged) bunching properties of bosons, i.e. the tendency of 
bosons to be found at the same position in space, and is defined as follows 

(2).Q^ ^ ^/vrf3x(v;t(x)v;t(x)^(x)v;(x)) 

^ [^^#x(^t(x)V;(x))]2 
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This quantity can be evaluated from the thermodynamic function F{Mq, Bq, V) 
by appropriate derivatives with respect to the chemical potential or with re- 
spect to the coupling strength g 0. 

At the critical temperature the dependence of this second order coher- 
ence factor on the scattering length a is depicted in Fig.§. It is apparent 
that with increasing scattering length g*^^^ decreases. This is consistent with 
the intuition that with increasing repulsive interaction bosons tend to avoid 
each other so that the probability to be found at the same position in space 
decreases. In the non-interacting case, i.e. for a = and fJ, = 0, this second 
order coherence factor assumes the value of 2 consistent with the well-known 
behaviour of a photon gas or chaotic field. 

Finally, let us consider the influence of confinement in one spatial di- 
rection on the second order coherence factor g^^-*. For this purpose we in- 
sert expression (EBf) for the density of states into the renormalization group 
equations (^21). The resulting dependence of g'-^-* on mean distance between 
the interacting bosons at fixed temperature is depicted in Fig.^ Tempera- 
ture and scattering length are chosen in such a way that Ath = 25a where 



-^th = y2T[fi /{mkBT) denotes the thermal de Broglie wave length. In the 
case of ^^Rb atoms, for example, with a scattering length of a = 5.3nm this 
condition corresponds to a temperature of T = 1.98^K. The full curve shows 
a case in which the characteristic length of confinement Lz is much larger 
than the thermal de-Broglie wave length. In this case the characteristic signa- 
ture of a second order phase transition is apparent at a critical scaled volume 
of mag nitude v = {V/N)X-^ = 0.456. As soon as the characteristic length of 
confinement Lz becomes comparable to the thermal de Broglie wave length 
this pronounced signature of a second order phase transition disappears. 



4 Summary 

The renormalization group constitutes a powerful method for evaluating ther- 
modynamic properties of interacting quantum gases beyond the limitations 
of mean-field theory. The previous considerations demonstrate that effects 
of confinement can be described by already existing momentum-space renor- 
malization group techniques as long as confinement can be modelled by pe- 
riodic boundary conditions. Within such an approach confinement leads to 
a modification of the density of the eliminated states in the Kadanoff trans- 
formation, thus affecting the flow equations of the renormalization group. 
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Typically, these modifications change the fixed points of the renormalization 
group equations and thus influence the critical properties in a signiflcant 
way. The presented treatment of effects of confinement by periodic bound- 
ary conditions constitutes a first step towards the more complicated final 
goal of obtaining an understanding of thermodynamic properties of trapped 
interacting quantum gases in realistic, smooth particle traps. 
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Figure 1: Flow of the scaled thermodynamic parameters {b{l),M{l)) for 
the ideal Bose gas according to the RG equations Eqs.(^. With increas- 
ing renormalization, i.e. increasing values of /, an initial point (60,^0) with 
60 = b{l = 0) ^ 1 and Mq = M{1 = 0) < in the thermodynamic parameter 
space is driven to higher effective temperatures, i.e. to smaller values of b{l), 
and to more negative effective chemical potentials. The unstable fixed point 
(6*, M^) = (0, 0) is indicated by a cross. 
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Figure 2: Schematic representation of the region in momentum space which 
is ehminated in the Kadanoff transformation. 
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Figure 3: Flow of the scaled thermodynamic parameters {M{l),G{l)) origi- 
nating in the RG equations Eqs.(0) for D = 3. The unstable, fixed point 
is indicated by a cross and the corresponding stable (s) and unstable (u) 
manifolds are indicated by arrows. 
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Figure 4: Dependence of the second order coherence factor g^'^'' on the scat- 
tering length a at the critical temperature. The mean distance between the 
bosons is denoted d = {V/Ny/^. 
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Figure 5: Isothermal dependence of the second order coherence factor (7^^^ 
on the scaled volume v = {V/N)X~^ for a thermal de Broglie wave length 
Ath = 25a. The continuous line corresponds to large confinement length 
Lz/\h ^ 1 and the dashed line to a small confinement length = 1.2Ath- 
For a large confinement length there is clear evidence of a second order phase 
transition at v = {V/N)X^^ = 0.456 which is smoothed out with decreasing 
confinement length. 
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